Constraining the evolutionary history of Newton's constant 
with gravitational wave observations 
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Space-borne gravitational wave detectors, such as the proposed Laser Interferometer Space An- 
tenna, are expected to observe black hole coalescences to high redshift and with large signal-to-noise 
ratios, rendering their gravitational waves ideal probes of fundamental physics. The promotion of 
Newton's constant to a time-function introduces modifications to the binary's binding energy and 
the gravitational wave luminosity, leading to corrections in the chirping frequency. Such corrections 
propagate into the response function and, given a gravitational wave observation, they allow for 
constraints on the first time-derivative of Newton's constant at the time of merger. We find that 
space-borne detectors could indeed place interesting constraints on this quantity as a function of 
sky position and redshift, providing a constraint map over the entire range of redshifts where binary 
black hole mergers are expected to occur. A LISA observation of an inspiral event with redshifted 
masses of 10 4 - 10 5 solar masses for three years should be able to measure G/G at the time of 
merger to better than 10 -11 yr _1 . 

PACS numbers: 04.30.-w,04.50.Kd,04.50.-h,04.80.Cc 
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I. INTRODUCTION 

The principle of equivalence is at the cornerstone of 
classical physics as evidenced by the opening paragraph 
of Newton's Principia Constructing a relativistic 

theory of gravity consistent with this principle was one 
of the key inspirations that led Einstein to general rel- 
ativity (GR) I^Tal- The Strong Equivalence Principle 
(SEP) states that all physics must satisfy: (i) the Weak 
Equivalence Principle (WEP), (ii) local Lorentz invari- 
ance and (iii) local position invariance. The WEP states 
that test-particle trajectories are independent of their in- 
ternal structure, while the latter two properties require 
experimental results to be independent of the 4- velocity 
and position of the laboratory frame. 

Many known alternative theories of gravity violate the 
SEP Q 1 , a classic example of which are scalar-tensor 
theories Q. These theories modify the Einstein- Hilbcrt 
Lagrangian by multiplying the Ricci scalar by a field <f> 
and adding a kinetic and a potential term (in the Jordan 
frame). Scalar-tensor theories are well- motivated from a 
theoretical standpoint, as they seem unavoidable in unifi- 
cation models inspired by particle-physics. In string the- 
ory, for example, there are numerous (dilatonic) scalar 
fields that couple to the scalar curvature, and perhaps 
also to the matter sector, then leading to additional vio- 
lations of the WEP [7Hll| . In loop quantum gravity, the 
Barbero-Immirzi parameter that controls the minimum 
eigenvalue of the area operator could be promoted to a 



field, leading to a classical coupling of Einstein gravity 
with a scalar- field stress-energy tensor [l2|, [l3| . 

From an observational standpoint, we also possess evi- 
dence for interesting physics that could point to a modi- 
fication of GR on cosmological scales, the prime example 
of which is the late-time acceleration of the universe. The 
cause of this acceleration is usually interpreted within the 
framework of pure GR as "dark energy" , a new form of 
matter that is or behaves similar to a cosmological con- 
stant. This form of matter has only recently (redshifts 
less than of order unity) become important in govern- 
ing the evolution of the universe. Another interpretation 
of the observed acceleration is that it is purely a grav- 
itational phenomenon, though the theory governing the 
dynamics of spacetime is not GR, but a new or modi- 
fied theory that has cosmological solutions that mimic 
a cosmological constant within pure GR. It is therefore 
of some interest to search for modifications to GR, in 
particular as a function of redshift. 

Often, gravitational theories that violate the SEP pre- 
dict a varying gravitational constant. In the scalar-tensor 
theory case, one can effectively think of the scalar as 
promoting the coupling between gravity and matter to 
a field-dependent quantity via G — > G(cj>), thus violating 
local position invariance as <p varies. Another interest- 
ing example are higher-dimensional, brane- world scenar- 
ios, where although physics remains Lorentz invariant in 
higher dimensions, "gravitational leakage" into the bulk 
inexorably leads to a time-varying effective 4D gravita- 
tional constant \A \. A constraint on the magnitude of 



1 The one known exception is Nordstrom's conformally-flat scalar 
theory @, which has been ruled out by light-deflection Solar 
System observations Q. 



this variation can also be interpreted as a constraint on 
the curvature radius of extra dimensions (l5| . 

Existing experiments searching for varying G can be 
classified into two groups: (a) those that search for 
the present or nearly present rate of variation (at red- 
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shifts close to zero); (b) those that search for secular 
variations over long time periods (at very large red- 
shifts). The most stringent bounds of the first class 
come from lunar ranging observations (l6j and pulsar 
timing observations [17|, [18 1 , but also include planetary 
radar-ranging and surface temperature observations 
of low-rcdshift millisecond pulsars [2(], [2l| . The strongest 
bounds of the second class come from Big Bang Nucle- 
osynthesis (BBN) calculations [22, [23| , but also include 
the evolution of the Sun [24[. The order of magnitude 
of the strongest constraint from either class is about 
G/G < 10~ 13 yr -1 , but it varies with the experiment. 

In this article, we show that gravitational wave (GW) 
observations of binary compact object coalescence can 
also provide constraints on the variation of G. We 
concentrate on these coalescences because of their high 
expected event rate and signal-to-noise ratio (SNR) 
with the proposed Laser Interferometer Space Antenna 
(LISA) [25j . We consider both nearly equal-mass, su- 
permassive black hole (SMBH) merger events (with total 
mass M ~ 1O 6 M0), as well as extreme-mass ratio, quasi- 
circular inspirals of a stellar-mass black hole (BH) into a 
SMBH. Only BHs are considered here because they lack 
internal structure that could lead to additional violations 
of the WEP [j| and modifications of the multipolar struc- 
ture of the gravitating body. From this standpoint, BH 
binaries are cleaner probes of fundamental physics. 

One interesting difference between previous tests of G 
and GW-based ones is that we are using black holes, 
ie. vacuum solutions, instead of neutron stars. As men- 
tioned earlier, scalar-tensor theories can be interpreted as 
variable- G theories, but in that case the variable G is re- 
ally only a variation in the coupling between gravity and 
matter. The gravitational constant G, however, serves 
more fundamentally as the quantity that defines the re- 
lationship between geometry (length) and energy, and 
such a quantity is not altered in scalar-tensor theories. 
We here consider the possibility that this fundamental 
constant varies with time. In this sense, we are looking 
at more generic variable- G theories, instead of restricting 
attention to theories that can reduce to a scalar-tensor 
theory. In Appendix [A] we discuss the relation between 
generic variable-G theories and scalar tensor ones in more 
detail. 

Another interesting difference between previous tests 
of G and the GW tests proposed here is that the latter 
constrain G/G at the time and location of the merger 
event. Therefore, once a sufficient number of events have 
been observed, a constraint map could be constructed, 
which would encode bounds on G/G along our past light- 
cone, ic. as a function of rcdshift z and sky position, 
out to a redshift of order 10 and limited mostly by the 
likelihood of merger events at that redshift. In contrast, 
the first class of tests discussed above can only bound 
G/G near redshift zero, while the BBN calculations only 
constrain linear changes in G from z = to z > 10 3 
and they arc degenerate with limits on the number of 
relativistic species. GW constraint maps would thus be 



excellent probes of local position invariance and the SEP, 
complementing previous bounds in a regime where the 
latter are inaccessible. 

Constraint maps are generic and applicable to any 
paramctrizablc GR deviation in the waveform, provided 
one possesses GW observations from sources at suffi- 
ciently large distances. The concept of GW tests, of 
course, is not new (see eg. [2rjW39| for a discussion of 
GW tests of Brans-Dicke and massive graviton theories), 
but until now constraint maps had not been really con- 
sidered. Such maps add a new dimension to the nature of 
GW tests; this is in contrast to most traditional tests that 
can only search for deviations in our immediate space- 
time neighborhood, or search for particular deviations 
that are by fiat constant over the age of the Universe. In 
this sense, GW constraint maps can provide invaluable 
global information on the SEP and local position invari- 
ance that is unattainable with Solar System observations. 

We find that a variability in Newton's constant intro- 
duces modifications to the phase of the GW response 
function of an inspiraling binary that could be measur- 
able by LISA to an interesting level. For comparable- 
mass BH inspirals with total redshifted mass M z = 
10 6 Af Q , one could constrain G/G < 10 _9 _yr _1 or bet- 
ter to redshift 10 (using an SNR of 10 3 Hcf), while for 
extreme mass-ratio inspirals with the same total mass, 
one could place roughly the same constraint out to red- 
shift 0.3 (using an SNR of 10 2 [H[). Such a bound be- 
comes stronger for intermediate mass BH inspirals, allow- 
ing constraints on the order of 10~ n yr^ 1 . Even though 
these bounds are, at a first glance, weaker than those 
from existing tests, the latter only constrain secular vari- 
ation in G, or oscillatory changes if we happen to be near 
a maximum of G today. GW observations arc in princi- 
ple sensitive to a much larger class of variations in G, 
if a sufficient number of events are detected to build a 
meaningful constraint map. 

We also find that there is an equivalence between the 
vacuum, varying-G model that we consider here and a 
time-varying chirp mass. Such an equivalence is estab- 
lished via the mapping (G/G)M — > M. In view of this, 
one can think of the tests presented here as probes of the 
time-constancy of the BH mass or horizons. Such time- 
variation could be induced, for example, due to gravita- 
tional leakage into the bulk in braneworld models [lj|. 
Therefore, the results presented here are truly a funda- 
mental GW test of the gravitational interaction. 

An outline of the remainder of the paper is as follows. 
Section [II] describes the parameterization chosen for the 
functional dependence Newton's constant. Section [TTT] de- 
rives how the inspiral is modified in post-Newtonian (PN) 
theory with a time- varying G to C(l/c 3 ), where c is the 
speed of light (this is commonly referred to as an expan- 
sion to 1.5PN order). Section HVl estimates the accuracy 
to which G/G could be measured by LISA, performing 
a simplified Fisher analysis and a time-domain dephas- 
ing study, using both frequency-domain PN waveforms 
and modified EOB, numerical templates. Section W\ con- 
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eludes with some discussion of future extensions of this 
work. In Appendix \K\ we further discuss the mapping 
between scalar-tensor theories and a variable-G theory, 
while in Appendix [B] we comment on possible degenera- 
cies due to cosmological expansion. 

Henceforth, we use the following conventions: we work 
exclusively in four spacetime dimensions with signature 
(—,+,+,+) Greek letters (/i, v>, . . .) in index lists 

range over all spacetime indices; Latin letters (a, b, . . .) 
in index lists range over all template parameter indices 
A a ; overhead dots stand for partial differentiation with 
respect to the time coordinate, E := dE/dt; the Einstein 
summation convention is employed unless otherwise spec- 
ified; finally, we use geometrized units G = c = 1 where 
possible, with c the speed of light in vacuum, but we 
shall mostly retain factors of Newton's constant G, as 
this quantity will be promoted to a time-function. 

II. PARAMETERIZATION OF THE 
VARIABILITY OF NEWTON'S CONSTANT 

In constructing the GW signal, we do not want 
to restrict attention to a particular alternative theory. 
Rather, we will assume that, during inspiral, the modifi- 
cation can be captured by simply promoting G to a func- 
tion of time in the GR equations of motion: G G(t). 
This is analogous to promoting G to a time-function 
in Newton's second law, a topic much studied in the 
past [43M46I ] . The G- varying framework we consider here 
can thus be thought of as an effective theory that could 
represent any alternative theory that leads to such a 
modification. In this section we discuss how we parame- 
terize G(t). 

At a given merger event, we will simply describe the 
evolution of the gravitational constant by 

G(t,x,y,z)nG c + G c {t c -t), (1) 

which is a Taylor expansion of G(t,x,y, z) to first or- 
der in time t and position (x, y, z) about the coales- 
cence event (t c , x c ,y c , z c ), where G c = G(t c ,x c ,y c , z c ), 
and G c = (dG/dt)(t c ,x Cl y c , z c ). Such an expansion is 
only valid provided the spatial variation in G is much 
smaller than the temporal variation, and the character- 
istic period of the temporal variation is longer than the 
observation window, the latter of which is on the order 
of a few years for LISA. 

Previous bounds on G(t) have assumed linear devia- 
tion from the present day value. In turn, this allowed 
BBN calculations to place such stringent constraints on 
G, by exploiting the huge "lever arm" provided by in- 
tegrating over cosmic time. However, any such bound 
will be largely insensitive to oscillatory components in G 
that have periods much less that the cosmic observation 
time. With the parameterization chosen in Eq. (fTJ), such 
oscillations can be captured, provided their frequency is 
much smaller than the inverse of the LISA observation 
time. 



Before proceeding, one might wonder how a variable 
G can lead to modifications to a vacuum spacetime, such 
as a BH inspiral. Newton's constant appears explic- 
itly on the right-hand of Einstein's equations, G = 
8irG/c 5 T^v, as a coupling constant between geometry, 
encoded by the Einstein tensor G tlVl and matter, given 
by the stress energy tensor T^ v . If T^ v = 0, as is the 
case in vacuum, then one might naively conclude that a 
variable G would not introduce any modification. In fact, 
this is essentially the case in scalar-tensor theories (see 
Appendix [A"]) . where the no-hair theorem is protected in 
vacuum by the scalar field's evolution equation. 

Newton's constant, however, plays a much more fun- 
damental role than merely a coupling constant: it de- 
fines the relationship between energy and length. For 
example, in the vacuum Schwarzschild solution, the re- 
lationship between the radius R of the black hole and 
the rest-mass energy E of the spacetime is R = 2GE/c 4 . 
Similarly, in a binary BH spacetime, the BHs have intro- 
duced an energy scale into the problem that is quantified 
by a specification of Newton's constant 2 . In our PN im- 
plementation of a varying G, these issues are avoided as G 
explicitly couples to a point-particle stress-energy tensor 
(the PN representation of a BH), and thus, it explicitly 
appears in the equations of motion. Henceforth, wc treat 
the variable-G modifications as induced by some effective 
theory and we do not concern ourselves with the specific 
model from which such modifications arise. 



III. POST-NEWTONIAN EXPANSION OF THE 
MODIFIED GW SIGNAL 

In the above parameterization of G(t,x,y, z) we have 
assumed that the spatio-temporal variation of G is suf- 
ficiently long to only induce a secular evolution over the 
duration of the observed inspiral and merger. Thus, the 
effect on the waveform will build up slowly over many cy- 
cles. We can then ignore the relatively short coalescence 
and ringdown phases of the merger to a good approxima- 
tion, and concentrate only in the inspiral. Modifications 
to the inspiral waveform can then be analyzed within PN 
theory; here we shall carry all expansions out to C(l/c 3 ) 
(ic. 1.5PN order). 

The time-dependence of Newton's constant modifies 
the generation of GWs due to corrections to the balance 
law during inspiral. Consider then a binary BH system, 
represented by PN point-particles with binding energy 
E = —rjMc 2 x/2, where rj = mim2/M 2 is the symmetric 
mass ratio, M = mi +7712 is the total mass, with mi ; 2 the 
binary's component masses, and x = (2nGMF/c 3 ) 2 ^ 3 , 



2 In a full numerical implementation, one could effectively model a 
time- varying G by introducing a time- varying length-scale, which 
is most easily done via a time- varying speed of light e(i). One 
must point out, however, that such modifications are unrelated 
to the variable speed of light theories of [47I , |4S(| . 
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with F(t) the orbital frequency. The balance law states 
that the rate of change of the binary's binding energy 
E is equal and opposite to the GW luminosity £ GW = 
32c 5 ?y 2 a; 5 /(5G), where the latter can be obtained from 
Isaacsson's stress-energy tensor [49|, The promotion 
of G — > G(t) in the conservative energy or Hamiltonian 
is comparable to a modification to Newton's second law, 
while this promotion in the GW flux is a modification to 
the dissipative sector of the theory. 

With the parameterization of Sec. HU the balance law 
can be rewritten as the following, modified differential 
equation for the orbital frequency: 



which then renders Eq. ([2]) into 
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where M. = rp/^M is the chirp mass. By construction, 
this expansion is valid to 1.5 PN order, ie. for orbital 
velocities v -C c to relative 0(v 3 ). This expression can 
be simplified by using the zeroth-order in G c relationship 
between time and frequency: 



t GR (F) = t c - — Gc Mx-S~ 8/5 [1 
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The promotion of G to a time-function with M con- 
stant is similar to the promotion of M to a time-function 
with G constant. One can establish this similarity by 
noting that the evolution of the GW luminosity function 
and binding energy is governed by the PN quantity x, 
which is always given in terms of the product G M. If 
one were to assume that M varied but not G, one would 
obtained the same type of evolution equation for the fre- 
quency [Eq. (J2J)], but with a slightly different prcfactor. 
One can then qualitatively think of a varying G as a 
varying M, or equivalently a varying chirp mass. 

The differential equation in Eq. ([2]) can be solved per- 
turbatively in G C M < 1 up to errors of 0(G 2 M 2 ): 



t(F) = t c 
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where evaluated at G(t) = G c . One might also 

expect corrections to quadrupole emission (such as dipo- 
lar or monopolar emission) due to a varying Newton's 
constant. Such non-quadrupolar emission could only be 
sourced by energy carried away by the G field. However, 
if we expect modified G-theories to still have well-defined, 
conserved energies, then the effective energy associated 
with G would scale as G 2 , which is sub- leading. If this 
behavior is generic, we can ignore such corrections, and 
we do so henceforth. 

Note that the PN equations are ill-behaved at the time 
of coalescence t c . Though we will stop the integration of 
the equations before this time, one might be concerned 



about the validity of Taylor expanding Newton's constant 
about t c as in Eq. (JT|). This may have been an issue if 
we had taken into account possible back-reaction effects 
of the orbital motion onto the local value of G, which 
we are not, as just discussed. Hence, to this order, it 
is valid to picture G(t) as an independent, cosmological 
field, whose functional behavior is not tied to the orbital 
motion of the binary. 

A modification to the frequency directly affects the 
GW phase evolution and thus the response function. The 
Fourier transform of this function in the stationary-phase 
approximation [5 11 ] becomes h(f) = A(f) exp(7^'), with 
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where u = TtG c M z f is a dimensionless parameter with 
/ the GW frequency, M z = (1 + z)M is the rcdshiftcd 
chirp mass, and </> c is a constant phase of coalescence. 
Q(l,(3) is some function of the inclination angle t and 
the polarization angle /3, which depends on the beam- 
pattern functions [51| • In Eq. © , we have omitted all the 
known PN amplitude corrections of GR and we have only 
explicitly written out the modifications to the amplitude 
induced by a time- varying G. 

The frequency dependence of the G c modification 
(/~ 13 / 3 ) is unique within the space of known alterna- 
tive theories of gravity [52| . Notice that the controlling 
factor of the expansion is changed by the G c terms from 
it" 5 / 3 to u -13 / 3 , which implies that to obtain corrections 
proportional to positive powers of the reduced frequency 
it, one would need to include corrections induced by 4PN 
terms. Comparing Eq. ([7]) to the analogous Brans-Dickc 
modification, we first notice that the G c term couples to 
M ?]~ 2 to leading order instead of to a positive power 
of the symmetric mass ratio 77, as in Brans-Dicke the- 
ory. In the latter, dipolar radiation activates when there 
are differences in the gravitational and inertial centers 
of mass (WEP violation), which requires objects with 
different self-gravitational binding energies (sensitivities) 
and mass ratios. On the other hand, the G c modification 
can be present in the absence of matter (ie. it modifies 
the point-particle representation of a BH spacetime). In 
spite of these differences, the leading-order G c modifica- 
tion can be recovered by the simplest form of the rece ntly 
proposed parameterized post-Einsteinian framework [52j , 
with the leading-order post-Einsteinian parameters 
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IV. PLACING A BOUND ON G 



With the insight gained from the analytic analysis of 
the previous section, we shall now concentrate on pro- 
ducing estimates of the accuracy to which G c /G c could 
be measured with a LISA GW observation. In Sec. IIV Al 
we briefly describe the basics of a dephasing accuracy es- 
timate and of a Fisher estimate. In Sec. IIVB1 we then 



proceed to calculate estimates using the PN expansion 
found in Sec. IIIII The resultant expressions, computed 
using only the leading order PN terms, give insight into 
the scaling behavior of the estimates with system param- 
eters. More accurate estimates can be obtained if one in- 
cludes higher order PN terms in the phasing. In Sec. IIV 01 
we compute the bounds on a G c /G c measurement using 
numerically generated waveforms that employ a certain 
resummation of the PN equations of motion and include 
higher order PN terms. We discuss all the different esti- 
mates in Sec. IIVDI 



A. Accuracy Estimates 

One criterion that can be used to estimate the accuracy 
to which G c can be measured is a dephasing measure. 
This measure requires that one compute the GW phase 
difference between a waveform with G c = Gq and G c = 
0. When this phase difference is equal to the inverse of 
the SNR 1/p, then one estimates that a G c correction 
of magnitude Go or larger could be detected. Such a 
measure follows from the analysis of [5^, [54[ . 

A more popular measure for determining the accuracy 
to which parameters could be estimated given a GW ob- 
servation is via a Fisher analysis. Consider then the 
Fourier transform of the response function h(f, A a ) as a 
frequency-series that depends on certain parameters A a 
that characterize the system. In our case, these param- 
eters define a 6-dimensional space, as we have neglected 
spin, that is spanned by A a = (Aq, t Cl 4> c , M, 77, G c ), 
where Aq is the first term of Eq. without the /~ 7//6 
factor. In reality, there are additional (so-called extrin- 
sic) parameters that also enter the response function and 
characterize the motion of the detector around the Sun 
through time-dependent beam-pattern function. One of 
the approximations we shall make here is to average over 
all angles so that the detector's proper motion can be 
ignored. 

In the Fisher scheme, the accuracy to which a template 
parameter A" can be measured is given by AA" = VE™, 
where E™ n is the (n, n) element of the covariance matrix 
(no summation over n implied) , defined as the inverse of 
the Fisher matrix 
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with the asterisk denoting complex conjugation. The 
quantity S n (f) is the spectral noise density, which de- 
scribes the noise in the observation. We shall here em- 
ploy the S n of (32I ] , which includes white-dwarf confusion 
noise. The upper limit of integration is chosen to 
be the minimum of the innermost stable circular orbit 
(ISCO) frequency (d^TrMg)' 1 and 1 Hz, while for the 
lower limit of integration /i Q we choose the maximum 
of the frequency 1 year prior to reaching the ISCO and 
10~ 5 Hz (see eg. [32| for further details). 

Another approximation we make here is to ignore cor- 
relations between G c and other parameters. With this 
approximation at hand, the Fisher matrix becomes diag- 
onal and its inverse tells us how well parameter A" can be 
measured: AA" ~ (r„„) ' . In our case, if no G c mod- 
ification is observed, we should be able to place a con- 



straint of approximately AG C 



< T, 



-1/2 



Gc/G c ^ v ± GcGc 

with the considcra- 



The estimate obtained from vE* 
tions described above, is what we refer to as a simplified 
Fisher analysis in this paper. 

Correlations generically tend to weaken the accuracy 
to which parameters can be measured. In our case, the 
inclusion of the mass ratio seems to deteriorate the con- 
straint by less than one order of magnitude. Other physi- 
cal effects, such as spin and eccentricity could also deteri- 
orate the constraint, but it is unclear whether this is the 
case. This is because such effects appear at high PN or- 
der, multiplied by a high power of the GW frequency and 
become important in the late inspiral, while the G c cor- 
rections modify the leading-order Newtonian term, multi- 
plied by a different power of the GW frequency, and dom- 
inate during the very early inspiral. Only a full Fisher 
analysis will reveal how much degeneracies deteriorate 
the constraint, but we relegate such considerations to fu- 
ture work. 



B. Stationary Phase, 1.5 PN Estimates 



where M z = M(l + z) is the redshifted total mass. In 
turn, this bound leads to the estimates 
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Clearly, the longer we observe the signal and the higher 
its SNR, the stronger the bound we can place, eg. a 3- 
year observation or an order of magnitude decrease in the 
total mass gains us roughly an order of magnitude on the 
constraint for each. 

Let us now perform a simplified Fisher analysis with 
the waveform of Sec. . To obtain less complicated 
expressions where the scaling with system parameters is 
apparent, we employ the so-called restricted PN approx- 
imation, where one neglects all amplitude corrections to 
the waveform. We find that the accuracy to which G c 
can be measured is approximately given by 
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where again we have retained only the leading-order G c 
correction to the GW phase in Eq. ([7]). In Eq. (fl~3|) . we 
have defined the reduced moment J p = I p /Ir and the 
moment of the distribution 1551 
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with the dimensionless variable x = f/fo and the ar- 
bitrary normalization constant /o = 10~ 4 Hz. Equa- 
tion (Ti"3"|) can be used to derive the following scaling rules: 



Let us now apply these measures to the frequency- 
domain PN waveform found in Sec. Mil We begin with 
the dephasing measure and make use of the fact that the 
GW frequency to leading PN order and zeroth-order in 
G c is given by 



c3/8 

/= ^—M^ 8 T- 3 ^ + 0(G c ), (10) 

57T 



where T is the integration time. Inserting this into 
Eq. l[7|). retaining only the leading-order G c term and 
requiring that *(G C ) - *(G C = 0) = 1/p leads to the 
bound 



(ag c ) 



< 



4 x 1(T 8 /100\ / i] 



EMRIs 



yr V 9 ) V 10 



(io ( 



10/3 



M 2 



(io ( 



10/3 



1-5 



( AG ) < 2 x ltr ' 5 ( * ^ 2 

V /smbh ~ yr \ 9 J \l/4 / 



(15) 



The scaling of this bound with total mass and mass ratio 
is different to that of Eq. (fT2"j) , because we have dropped 
in Eq. (fT5"j) the J33 dependence, which itself depends on 
observation time, total mass and mass ratio. 

It is somewhat surprising that the Fisher and dephas- 
ing estimates are quite constistent for EMRIs, though 



7 



much less so for SMBH mergers. One reason for this is 
that we have neglected sub-leading PN corrections in the 
frequency-domain waveform of Eq. 0, which one would 
expect to have a much larger effect for equal mass sys- 
tems. To check this, we can compute the sub-leading PN 
corrections to Tq q , which leads to the bound 



AG, < 



65536 1 V 2 13/3 1 

25 p M z V ^33 

743 22 \ 2/3 J 31 

126 -YV ifoM7:) J33" 

O ^33 



(16) 



For a system with the parameters M z = 10 6 M Q , r/ = 1/4 
and p = 1000, Eq. (pH) gives AG C < 8x 10~ 7 yr" 1 , which 
is more than an order of magnitude stronger than that 
obtained in Eq. (JTSJ). This result gives us a sense of the 
importance of the inclusion of higher-order corrections to 
the GW phase and it indicates that the SMBH bound in 
Eq. (|15p may be unreliable (note that due to the rj depen- 
dence in (|16p the EMRI estimate is largely unaffected). 
With this in mind, the dephasing and simplified Fisher 
measures suggest that roughly the same order of mag- 
nitude constraint on G c from either EMRIs or SMBHs 
could be obtained by LISA. 



C. Effective-One-Body Estimates 

In the previous subsection we have seen that the in- 
clusion of higher-order corrections can influence the es- 
timated accuracy to which G c can be measured. For 
this reason, in this subsection we calculate a more accu- 
rate GW phase, obtained via numerical methods through 
the so-called effective-one-body (EOB) formalism [5g- 
|69| . Essentially, one numerically evolves the bodies' tra- 
jectories by solving the Hamilton- Jacobi equations for 
a certain resummed, high-order PN Hamiltonian. The 
waveform is then constructed via another resummed PN 
expression that depends on the evolved trajectories. We 
shall not here review the EOB framework, but instead 
we follow the prescription described in detail in (70j . In 
particular, we employ the uncalibrated 5.5PN order EOB 
model, including ^-deformations of [7(1 IzJ > although we 
find our results are rather insensitive to the specific EOB 
model chosen. 

The G c modification is straightforward to implement 
in the EOB scheme. The conservative dynamics (the 
Hamiltonian) depend on reduced coordinates and con- 
jugate momenta (denoted by overhead hats), which are 
normalized to the total mass and mass ratio p = rjM re- 
spectively. The G c corrections to the Hamiltonian can 
then be modeled by rescaling certain reduced quanti- 



ties by G(t) (see eg. Sec. Ill of [JJ) 3 : t -> i MG{t), 
r -> f MG(t), p,/, -> MpG(t). The dissipative dy- 
namics (the radiation-reaction force) can be modeled by a 
certain Pade resummation of the energy flux, which itself 
is a function of the orbital velocity. The G c corrections 
to the radiation-reaction force can then be modeled by 
reinstating the factors of G and employing Kepler's third 
law, which states v = [G(t)Mw] 1 / 3 , ie. the ^-component 
of this force is rescaled via T<j,{v) — > J r 0[G(t) 1 / 2 u] (see 
eg. Sec. Illb in In this way, the Hamiltonian and 

the dissipative force become explicit functions of time, 
leading to modified Hamilton- Jacobi equations with new 
terms proportional to G c , whose numerical solution al- 
lows us to compute G c -modificd EOB waveforms. 

We begin first by applying the dephasing measure to 
these EOB waveforms. Figurc[T]plots the difference in the 
GW phase as a function of time in units of months for dif- 
ferent values of G c . Clearly, as G c is chosen to be larger, 
the difference in GW phase also increases. With the 
criteria described above, these evolutions suggest LISA 
would be able to detect a G c ~ 10 -9 yr _1 when observ- 
ing extreme-mass ratio inspirals (EMRIs) with p ~ 250 or 
SMBH mergers with p ~ 3300. As expected from Eq. (J7J), 
EMRI systems are preferred over SMBH/SMBH coales- 
cences due to the extra factors of rj multiplying the G c 
correction (or intuitively, due to the slower inspiral rate, 
an EMRI spends many more cycles in the strong-field 
regime prior to merger, enhancing the secular deviations 
induced by G c ). 

Let us now apply a numerical version of the simplified 
Fisher analysis to the modified EOB waveforms. Via Par- 
seval's theorem, we can convert a frequency-domain inte- 
gration into a time-domain integration. The Fisher ma- 
trix, however, contains the spectral noise density in the 
denominator, which makes this conversion non-trivial. 
We shall here treat this density as an averaged quantity 
that we evaluate at the time t(f). With this in mind, 
we can rewrite the (G c , G c ) element of the p 2 -normalized 
Fisher matrix as 



F 



and the SNR is simply 



lyr 



dt 



Sr. [/(*)] 



dh 



dG c 



lyr 



dt 



Sn[f{t) 



\h\' 



(17) 



(18) 



The pseudo-averaging of the spectral noise density is sim- 
ilar to what is done in the stationary phase approxi- 
mation, where the generalized Fourier integral that de- 
fines the Fourier transform is evaluated at the stationary 



3 Since we study quasi-circular orbits, we can neglect the 8 evo- 
lution. The Hamilton-Jacobi equations do not acquire overall 
multiplicative G c corrections, as this can be cast completely in 
terms of dimensionless quantities. The only G c corrections arise 
via the rescaling of the coordinates and momenta inside the ef- 
fective Hamiltonian. 
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FIG. 1: Difference in the GW phase between a G c = evoiution and a G t / evolution, as a function of time in units of 
months. Left panel: fiducial extreme-mass ratio inspiral with r\ = 10~ J and redshifted mass M z := M(l + z) = 10 6 in the last 
stages of coalescence (the orbital radius r/M £ (6, 10)). Right: fiducial SMBH inspiral 77 = 1/4 and M z — 10 6 in the last stages 
of coalescence (the orbital radius r/M £ (6,95)). 



point 4 . The accuracy to which G c can be measured is 
then simply AG C < p-^f^J" 1 / 2 . 

We compute the relevant element of the Fisher ma- 
trix numerically. The partial derivatives in Eq. (|17p are 
approximated via the rule dh/dG c ~ [h(Go) — h(Q)]/Go 
for some very small Go- We have calculated the Fisher 
estimate for different values of Go € {10 -10 , 10 _7 }yr _1 
and found the constraint to be robust to such changes. 
The waveforms employed are aligned in an initial time 
window, in the low-frequency domain: we have time and 
phase-shifted them so that they agree to at least one part 
in 10 6 during the first day of evolution. Such a calibration 
amounts to a matched filtering maximization of the SNR 
in white noise over the time and phase of coalescence [§9[ 
(see [7l[ for more details on this procedure). 

Using this numerical Fisher prescription and the mod- 
ified EOB code, we find that for an equal-mass, SMBH 
inspiral with M z = 10 6 A/ Q , G c could be measured 
to AG C < 3.5 x 10- 9 yr- 1 (1000/ / o). Similarly, for an 
extreme-mass ratio inspiral with the same total mass 
and mass ratio 77 = 10 -5 , G c could be measured to 
AG C < 4.5 x Hr 9 yr -1 (100/p). These estimates are in 
rough agreement with the graphical dephasing measure 
presented in Sec. IIV C[ as well as with the frequency- 
domain simplified Fisher analysis of Sec. IIVB1 keeping 
in mind the various approximations employed. 



4 Ref. [73(| has verified that this time-domain implementation of 
a Fisher analysis produces very similar results to the standard 
frequency-domain implementation, when the Fourier transform 
of the waveform in the latter is treated in a stationary-phase 
approximation. 



D. Comparing Estimates 

We have found that the dephasing and simplified 
Fisher measures are rather consistent among themselves. 
Differences, however, arise when considering a Taylor- 
expanded, stationary-phase approximated, 1.5PN order 
waveform versus a modified EOB wave. These similari- 
ties and differences can be better appreciated in Table Q] 
At some level, we expect the EOB waveforms to lead to 





PN deph. 


PN Fisher 


EOB deph. 


EOB Fisher 


EMRIs 
SMBH 


3 x 10" s 
1.5 x 10" 7 


4 x 10" s 
8 x 10" 7 


2.5 x 10" a 
3.3 x 1(T 9 


4.5 x 10" M 
3.5 x 1(T 9 



TABLE I: Constraints on G c /G in units of yr -1 , for a 1 year 
observation, using the the dephasing (first and third columns) 
or the Fisher measure (second and fourth columns). The first 
row assumes an EMRI with total redshifted mass M z = 10 6 , 
symmetric mass ratio 10 -5 and SNR of 100, while the second 
row assumes an equal-mass SMBH inspiral with the same to- 
tal redshifted mass and SNR of 1000. In the case of the PN 
Fisher estimate for SMBH inspirals, we have used the higher- 
order expressions, as the leading-order term seems unreliable. 

a better bound, as they include many more G c correc- 
tions than the analytic 1.5PN model. As already found 
in Sec. IIV B[ these higher-order PN terms have more of 
an impact on comparable mass systems than on EMRIs. 

Overall then, the numerical estimates suggest a typi- 
cal GW inspiral could be used to constrain G c /G c to a 
level of ~ 10~ 9 yr -1 . The best bounds come from in- 
termediate mass BH inspirals, eg. for an inspiral with 
redshifted masses (10 4 : 10 5 )M Q , one can constrain 
G c /G c < 10- 11 vr -1 , assuming an SNR of 10 3 and an 
observation time of 3 years. Such a constraint may even 
become stronger once spin and eccentricity are included, 
as then one can possibly obtain more cycles in the sen- 
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sitivity band. Realistically, however, although the most 
optimistic systems begin to become comparable to So- 
lar System constraints, typical SMBH inspirals at high- 
redshift will probably not lead to such strong constraints. 
Nevertheless, as explained in the Introduction, a suffi- 
ciently large number of observations will allow for the 
construction of a map of the constraints on G as a func- 
tion of redshift and sky position, a task that is impossible 
with Solar System observations. As such, GW tests con- 
stitute a new, complementary class of G tests, allowing 
for a unique opportunity to build a constraint map during 
the cosmological epoch of structure formation. 



V. CONCLUSIONS 

The promotion of Newton's constant to a time- 
dependent quantity introduces modifications in the GW 
frequency and phase of an inspiraling binary. We have 
analytically computed the G modifications to the re- 
sponse function in the stationary phase approximation 
to 1.5PN order. We have also numerically solved a G- 
modificd version of the EOB equations of motion to es- 
timate the magnitude of the G corrections to higher PN 
order. We have used both a dephasing and a Fisher mea- 
sure to obtain order of magnitude estimates of the pos- 
sible constraints that could be placed with LISA given a 
GW observation. 

An observation of one merger event could be used 
to constrain deviations from G = const, at the space- 
time location of merger. A single GW observation of an 
equal-mass inspiral event with total redshifted mass of 

10 5 M Q for 3 years should be able to constrain G c /G c to 
10 -n yr -i 

Given a sufficient number of GW observations from 
events at different sky locations and redshifts, one could 
in principle construct a constraint map of G as a function 
of sky angles and redshift. Conversely, if G is in truth 
not constant, and deviations are measurable, GW obser- 
vations could allow for a reconstruction of G(t,x,y, z) 
along our past lightcone. Solar system bounds constrain 
G at redshift zero, while cosmological BBN calculations 
bound a secular drift in G from redshift to z > 10 3 . 
With these restrictions, existing bounds are a couple of 
orders of magnitude stronger than putative GW obser- 
vations; however, GWs will allow for entirely new con- 
straints on G at intermediate redshifts, and thus as more 
events are observed throughout the universe observations 
will become increasingly sensitive to oscillatory compo- 
nents in G (with periods greater than of order a year). 

Certainly, the same techniques described here could 
be used to place constraints on G from ground-based 
GW observations of solar mass BH mergers with 
LIGO/VIRGO/GEO. However, the lower expected SNR 
of these events, together with the shorter integration 
time, imply the bounds would be much weaker than ex- 
isting bounds or putative LISA bounds. Furthermore, 
since the typical, observable solar mass merger events 



will occur at relatively low redshifts, ground-based de- 
tectors will also not be able to offer a measurement that 
is complimentary to the already strong solar system con- 
straints. 

Future work could concentrate on performing a more 
detailed data analysis of the accuracy to which G c can be 
measured, in particular taking into account possible de- 
generacies with other GR binary merger parameters. We 
have estimated that such correlations will deteriorate the 
accuracy to which G can be constrained by roughly an or- 
der of magnitude, when degeneracies with the chirp mass 
and symmetric mass ratio are included, a result in qual- 
itative agreement with j32|. The inclusion of additional 
physical effects, such as spin, eccentricity or higher-order 
PN harmonics, however, might break some of these de- 
generacies and restore some of the strength in the con- 
straint [3714391 ] . The spin corrections, for example, first 
enter in the phase at 1.5 PN order via a spin-orbit interac- 
tion, and as such are only dominant close to merger, while 
the G modification also affects the early inspiral. Even 
if these degeneracies deteriorate the bound dramatically, 
however, the GW constraint remains unique in its sensi- 
tivity to non-linear behavior in G, if a sufficient number 
of detections are made to measure G at different locations 
in the universe. 

Another interesting avenue to pursue would be to 
investigate how the bounds on a varying G could be 
strengthened if in addition to the GW signal an electro- 
magnetic (EM) counterpart was observed. When deal- 
ing with pure vacuum spacetimes, such an EM signature 
is not likely, but recent work on numerical simulations 
of mergers embedded in circumbinary disks suggest that 
EM signals might be emitted [741 - 1791 ] . Given a coincident 
EM/GW detection, one would then be able to determine 
the position of the source [8(| and its redshift (8ll - [83l |. 
while given a cosmology, in the form of a distance- redshift 
relationship, one would then be able to independently 
estimate G at the source. The tandem GW/EM obser- 
vation would thus simultaneously constraint G and G at 
the merger event. 

One last path to explore is the possible degeneracies 
between the G c effect and the GR time-dependence of 
redshift z, induced by the expansion of the Universe. 
Both effects lead to a modification to the GW phase that 
is proportional to the same power of the GW frequency, 
making these effects exactly degenerate. We performed 
a preliminary exploration of such a degeneracy in Ap- 
pendix [B] where we found that when cosmological infor- 
mation is employed, the z effect can be subtracted off to 
still allow for bounds on G at the levels discussed in this 
paper. 

Note added after submission: Shortly after this 
paper was submitted, one [&H appeared in preprint con- 
sidering the possibility of constraining the size of string- 
theory inspired extra-dimensions, which induce a time- 
dependence in the BH mass. As we discussed earlier, 
such an effect is equivalent to a varying G, with the map- 
ping G c /G — > M/M. Because the extra-dimension effect 
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is strongest for the smallest BHs, we note that EMRIs 
are the best systems to explore this effect with GWs. A 
direct GW detection could place a bound on the size of 
the extra-dimension of I > 1 micrometers, assuming a 
constraint on M/M on the order of 10 -9 yr -1 . Such a 
constraint can be improved by one order of magnitude if 
the system is observed for more than one year. 
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where r — \J x 2 + y 2 + z 2 . From the above, the identifi- 
cation of an effective G via the action gives 



G B D 



G 



2GM r- 1 (3 + 2w E 



(A4) 



Clearly, G bd varies with distance from the source and 
it asymptotes to G BD ~ G (f)^ 1 in the limit M/r <C 1. 
However, from the (0, 0)-component of the metric, the 
effective Newtonian potential is — G G ffM/r, with 



Gog 



G r 



1 + (3 + 2lu bd ) 



(A5) 



Appendix A: Generic Variable-G Theories and 
Scalar- Tensor Models 



In this appendix, we discuss the relation between the 
generic variable- G theories we consider in this pa per and 
one particular scalar-tensor theory: Brans-Dicke [851 ] . In 
Brans-Dicke theory the action is given by 



S = 



16ttG 



Ft 



Sri 



(Al) 



where g is the determinant of the metric, R is the Ricci 
scalar, is a scalar field, w BD is a coupling constant, 
commas stand for partial differentiation and S m &t. stands 
for additional matter degrees of freedom. This action 
leads to modified field equations, as well as an evolution 
equation for the scalar field: 



9V 3 + 2w BD V 30 



(A2) 



where O g is the D'Alcmbertian operator of curved space- 
time and T is the trace of the matter stress-energy tensor. 

Brans-Dicke theory is often interpreted as a variable-G 
theory by defining a new Brans-Dicke gravitational pa- 
rameter via G BD = G _1 . Thus, if varies in spacctime, 
so does the effective coupling parameter G BD . However, 
note that Newton's constant G is still exactly that: a con- 
stant. What is changing is the coupling between other 
forms of matter and gravity, as mediated by 0. 

In order to illustrate this point, consider the original 
weak-field solution found by Brans and Dickc describing 



We notice that this effective potential is the same as 
that of general relativity, but with a constant rescaling 
of G — > G ff. Note further that, in the spatial com- 
ponents of the metric, a different constant rescaling of 
G would bring those components to a form analagous 
to Schwarzschild. One of the consequences of this, as 
orginally pointed out by Brans and Dicke, is that to 
this order in the weak-field approximation circular orbits 
are completely unaffected compared to GR predictions in 
the solar system, although light deflection is for example 
modified. 

The example above shows that Brans-Dicke theory is 
not fundamentally a variable G theory: absorbing into 
an effective variable-G coupling at the level of the action 
is a formal device and care must be taken in how this is 
interpreted. More significantly, Brans-Dicke theory docs 
not predict any GR deviations in vacuum. In particu- 
lar, binary BH mergers, as considered in this paper, are 
completely unaffected. The reason for this is that Brans- 
Dicke theory (and all similar scalar-tensor theories) only 
alter the coupling between matter stress-energy and grav- 
ity; such theories do not change how the stress-energy of 
the gravitational field itself back-reacts on the geometry, 
which is still controlled by (the unsealed) G. 

We can elucidate this further by taking an effective 
field theory point of view. Consider the GR action, lin- 
earized in the perturbations h^ v about Minkowski space- 
time 77^, ie. g^j = r\^ v + h^: 



S = 



d 4 x [(h a " tap - r) ap O v h ap ) 

+ 167rG ma tt£matt 

+ a- 1 (d^Gg) (d»G g 



167tG s £ gw 



(A6) 
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where is the D'Alcmbertian operator of flat space, 
£matt is some matter Lagrangian density and £ GW is 
the second-order, effective Lagrangian density associated 
with GWs. We have included two coupling constants 
here: G ma tt associated with the coupling between gravity 
and matter; and G g associated with the self-gravitational 
coupling. In GR, both of these coupling constants are 
identical, but in Brans-Dicke G ma tt = G/<j>, whereas G g 
is still G. Essentially, in Sec. IIIII of this paper we have 
promoted G g to a function of spacetime, and this can- 
not be mimicked by Brans-Dicke or similar scalar-tensor 
theories. 

The last term in Eq. (|A6|) gives dynamics to the grav- 
itational constant, which leads to the equation of motion 



(A7) 



where a is a coupling constant. Notice that this cou- 
pling is different from that which arises in scalar-tensor 
theories, as in the latter the scalar field couples to the 
trace of the stress-energy tensor, instead of to GW non- 
linearities. 



Appendix B: Cosmological Degeneracies 

The modification to the GW response function, de- 
rived in Eq. ([Jj), is not unique. As mentioned in the 
Introduction, this effect can be also obtained in a the- 
ory where G is not varying, but Ai is time-dependent. 
Both of these are wave generation effects, where the evo- 
lution equation for the GW frequency is modified due 
to the new time dependence. In addition to this, a sim- 
ilar modification arises if the background cosmology is 
time-dependent, ie. in an expanding universe. Although 
such a modification is a wave propagation one, it leads to 
the same frequency modification in the frequency-domain 
GW phase to leading order in the stationary-phase ap- 
proximation. In this appendix we show how this comes 
about and consider its implications in the constraint of 
G c . 

Consider an expanding cosmology, where the metric 
upon which GWs propagate, is described by an unper- 
turbed Fricdmann- Robertson- Walker line element: 



ds 2 



-dt 2 + a 2 \t)dx l dx 3 Si 



(Bl) 



where a(t) is the scale factor and t is cosmic time. The 
scale factor allows us to define the Hubble parameter 
in the usual way: H(t) = a/a, where the overhead dot 
stands for partial differentiation with respect to cosmic 
time. Using the standard definition for the redshift, 
z = ao/cL—1, where do = a(to) corresponds to the present 
value of the scale factor, we can derive the evolution equa- 
tion 

Ay 

- = (l + z) 2 H -(l + z)H(z), (B2) 

The constant H is the value of the Hubble parameter 
H(z = 0) today, H = (70.5 ± 1.3) km s" 1 Mpc -1 



Notice that this expression differs from that derived by 
McVittie HU [dz/dt = (1 + z) H -H(z)], because here 
we differentiate redshift with respect to cosmic time at 
the source, instead of at the observer. 

The Hubble parameter depends on the Fricdmann 
equations, and neglecting any G(t) dependence, one can 
write it as 



(B3) 



where Q m ,o and f2 c stand for the present-day matter den- 
sity and the cosmological constant respectively, while the 
. . . represent other terms that are negligible in the regime 
z < 10. In a matter-dominated era, the regime of most 
relevance for LISA, we can approximate the Hubble pa- 
rameter as H{z) ~ Hq (1 + z) 3 ^ 2 0^{q, thus leading to 



dz 
~dt 



Hn 



5/2 n l/2 



(i + z y -(l + zY'^n 



m.O 



(B4) 



Notice that this effect not only affects gravitational wave 
observations, but also any electromagnetic one that in- 
volves a sufficiently long, coherent time-integration. 

Let us now connect z to the GW response function. As 
already pointed out, the response function, as measured 
on Earth, is sensitive to the redshifted chirp mass M. z = 
(1 + z)Ai. Though since the redshift changes with time, 
which to leading order we can expand about the time of 
coalescence as 



z{t) 



z c + 



dz 
dt 



(t - t c ) : 



(B5) 



+ h q (i + Zc y - (i + Zc f 2 nlil (t-t c ) 



one sees that M. z becomes time-dependent. Such an ef- 
fect contributes to any cosmological observable that de- 
pends on redshift, as first pointed out by McVittie (88j . 

The calculation of the GW observable is not yet com- 
plete as one must relate the time dependence (t — t c ) in 
Eq. (|F35|) to GW frequency. As pointed out in Sec. IIIII 
the Fourier transform of the GW response function in 
the stationary phase approximation is given by Eqs. ©- 
J7J), except that now redshift has become a function of 
time in the inspiral frame. One can rewrite this function 
of time as a function of GW frequency, through Eq. §5§ 
with the substitutions F —> f/2 and M — >• M(l + z c ), 
which appropriately rcdshifts all physical lengths into the 
observer frame. 

Concentrating on the dominant term in the GW phase, 
and defining \&i. , to be the leading order term of Eq. ([7]) 
and S^z to be the new contribution to Eq. ([7]) , such that 
* = *i.o.(l + S^i), we find 



25 



(B6) 



where now u = nG c M Zc f and M Za = M(l + z c ). The 
new relative contribution, 8^ z , depends on /~ 8//3 , or 
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cquivalcntly, « -8 / 3 , just as the G c term in Eq. (O, and 
thus the two effects are degenerate. 

This raises the question of whether the z effect could 
overwhelm the G c effect, rendering the bounds derived 
in this paper unattainable. To answer this question, we 
will use the results of the PN dephasing/Fisher analysis 
of the G c effect in Sec. IIVBI Using the measured val- 
ues of the cosmological parameters obtained in the five- 
year WMAP analysis (13, H (H = 70.5 km s" 1 Mpc" 1 , 
^m,o — 0.26), and defining the ratio r = 5^> q of the 
correction induced by the i term to the 5^!^ correc- 
tion induced by the G c term (the first term in the second 
line of Eq. ©), we find that r « 0.008 (10~ 8 yr^/Gc) 
at z c = 1, growing to r w 0.02 (10~ 8 yr _1 /G c ) at z c = 5 
for equal-mass binaries. Since this ratio is less than one, 
it implies the cosmological effect is not significant, and 
does not need to be accounted for. 

One might wonder though why we did not use the more 
sensitive bounds for G c obtained from the EOB analysis 
in the above. The reason is that the z effect is purely 



a propagation effect, and would not be enhanced by in- 
clusion of higher-order PN terms. This is in contrast 
to the effect of G c , which changes the local dynamics 
of the binary, and for near-equal mass systems this is 
strongly reflected in the higher order PN terms. Nev- 
ertheless, even if we did assume the EOB G c bounds 
carried over verbatum to i, one can still subtract the 
z effect using cosmological data prior to testing for G c . 
How much this will reduce the z effect depends on how 
accurately the cosmology is known. Assuming a rela- 
tive error of 10% in Hq and £l m .o (consistant with how 
accurately these numbers are known today), the ratio 
is reduced to r « 0.2 (10~ 10 yr _1 /G c ) at z = 1 and 
r w 0.5 (10~ 10 yr _1 /G c ) at z = 5. From new and antici- 
pated observations, provided for example by Planck [ooj ] . 
it is not unreasonable to assume by the anticipated time 
of LISA the relative error in the cosmology would be an 
order of magnitude less, further reducing r by an order of 
magnitude. This is then well below our most optimistic 
bounds for measuring G c . 
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